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Elliptic flow and HBT are studied for non-central relativistic nuclear collisions. Azimuthal asym- 
metries show up in both elliptic flow and HBT radii and are calculated in both collisionless and 
hydrodynamic limits relevant for peripheral and central collisions respectively. Determining the 
reaction plane and measuring the HBT radii as function of the angle between the reaction plane 
and the particle momenta can determine the physical quantities as source sizes, deformations, emis- 
sion times, duration of emission and opacities. Comparison to SPS data and predictions for RHIC 
and LHC energies are given. The centrality dependence with and without a phase transition to a 
quark-gluon plasma is discussed - in particular, how the physical quantities are expected to display 
a qualitative different behavior in case of a phase transition. 
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I. INTRODUCTION 



The azimuthal asymmetries in momentum spectra usu- 
ally referred to as directed and elliptic flow are large in 
nuclear collisions at intermediate energies 0] which al- 
lows experimental reconstruction of the reaction plane. 
Also in relativistic nuclear collisions at the AGS Q and 
SPS H anisotropic flow has been found in semi-central 
collisions. Likewise, HBT analyses has successfully been 
made, however, so far without simultaneous reconstruc- 
tion of the reaction plane. The statistics may permit 
this in the near future and the prospects looks good for 
the RHIC and LHC colliders, where the multiplicities are 
much higher and the number of pairs entering HBT anal- 
yses grow with the multiplicity squared. Particle inter- 
ferometry (HBT |Q) may also show source asymmetries 
^,[1 spatially and complement the momentum space in- 
formation of flow. Earlier studies of directed and elliptic 
flow by hydrodynamics j7| as well as by the cascade code 
RQMD has been compared to AGS and SPS data 
with some success. The data shows interesting transverse 
momentum and rapidity dependences for pions and pro- 
tons that cannot be explained by these models. 

The aim of this work is to return to the basic physics 
that leads to azimuthal anisotropics in semi-central col- 
lisions in both coordinate and momentum space which 
are relevant for the HBT radii and anisotropic flow re- 
spectively. We shall derive simple analytical formulas for 
elliptic flow and HBT radii in both the collisionless and 
hydrodynamic limits relevant for peripheral and near- 
central collisions respectively. A combined analysis of 
HBT and anisotropic flow yields detailed information on 
the particle source at freeze-out as well as the expansion 



from the initial collision and up to freeze-out. General 
properties as the initial geometry, shadowing effects, ex- 
pansion and freeze-out can be determined through de- 
tailed measurements on source sizes, deformations, life- 
times, duration of emission, opacities, transverse and el- 
liptic flow. We discuss the effects of a phase transition 
occurring when the energy densities become high enough 
to form a quark-gluon plasma at some semi-centrality. 
The centrality, Et or dN / dy dependence of a number of 
measurable quantities as the ratio of the elliptic flow to 
the deformation, duration of emission, opacity, etc. may 
display an interesting behavior. 

The paper is organized as follows. First we describe 
the geometry of semicentral relativistic nuclear collisions 
and discuss the azimuthal asymmetries and deformations 
of the overlap zone in sec. II. In sec. Ill we calculate el- 
liptic flow in the collisionless and hydrodynamic limits 
relevant for peripheral and near central collisions respec- 
tively and compare to recent SPS data || . HBT radii are 
calculated for deformed sources in sec. IV. In sec. V wc 
discuss the qualitative behavior of a number of physical 
quantities from elliptic flow to HBT radii - with and with- 
out a phase transition to a quark-gluon plasma. Finally 
a summary is given. 



II. GEOMETRY OF SEMI-CENTRAL 
COLLISIONS 



In non-central nuclear collisions the reaction plane 
breaks azimuthal symmetry. The asymmetry decreases 
with centrality and vanish for the very central collisions 
(see Fig. 1). 
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We will use the term centrality as impact parameter b 
in the collision. It is not a directly measurable quan- 
tity but is closely correlated to the transverse energy 
produced Et, the measured energy in the zero degree 
calorimeter and the total particle rapidity density dN/ dy. 
The latter is again approximately proportional to the 
number of participating nucleons 
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The number of participants for colliding two spherical nu- 
clei of radius R is shown in Fig. 2 as function of impact 
parameter. 

The standard geometry has the z-direction along the 
longitudinal or beam axis and the x-direction along the 
impact parameter b. Thus (x,z) constitutes the reaction 
plane and (x, y) the transverse direction with y perpen- 
dicular to the reaction plane. It is convenient — also 
for comparison to HBT analyses — to employ gaussian 
parametrizations for anisotropic sources in both trans- 
verse directions 
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Here, the gaussian transverse sizes R x and R y increase 
with time as the source expands after the collision. The 
initial transverse radii, i.e., the rms radii of the over- 
lap zone weighted with nuclear thickness functions, are 
shown in Fig. 2. For HBT the relevant radii are, however, 
those at freeze-out which are larger due to expansion as 
discussed below. 

The azimuthal asymmetry or "deformation" of the 
source can be defined as the relative difference between 
the gaussian radii squared 
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A simple estimate can be obtained from the full trans- 
verse extent of the initial overlap of two nuclei of mass 
number A and radius Ra — 1.2A 1 / 3 colliding with impact 
parameter b (see Fig. 1). They are: R x = Ra — b/2 and 
Ry = yRjs, — b 2 /4, and the corresponding deformation 



is 



S = 



2R A 
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The rms radii of the nuclear overlap zone weighted with 
longitudinal thicknesses results in deformations (see Fig. 
2) that are slightly smaller than Eq. ([IJ) at semicentral 
collisions. However, as the source expands the deforma- 
tion 5 ~ (Ry — R x )/(R y + Rx) decreases for two reasons. 

Firstly, the expansion increase (R x + R y ) as is found 
in HBT and Coulomb analyses |11|. Secondly, (Ry — R x ) 
decrease because the average velocities are larger in the 
x- than y-direction. The latter is a consequence of the 



experimentally measured positive elliptic flow (v% > in 
Eq. (|J)) in relativistic nuclear collisions where shadowing 
is minor . Measuring the decrease of the deformation 
with centrality will reveal important information on the 
expansion up to freeze-out. For very peripheral collisions, 
where only a single nucleon-nucleon collision occurs, the 
source must be azimuthally symmetric, i.e., the deforma- 
tion must vanish and therefore Eq. (4), which assumes 
continuous densities, breaks down. 



III. ELLIPTIC FLOW 



The reaction plane, which breaks azimuthal symmetry, 
has been successfully determined in non-central heavy ion 
collisions from intermediate up to ultrarelativistic ener- 
gies H . The particle spectra are expanded in harmonics 
of the azimuthal angle <f> event-by-event H 
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where n(x,p) is the particle distribution function in space 
and time x — r,t), and 4>r is the azimuthal angle of the 
reaction plane. Assuming that the experimental uncer- 
tainties in event plane reconstruction can be corrected 
for, each event can be rotated such that 0# = 0. The 
asymmetry decreases with centrality (see Fig. 1+2 and 
Eq. (0)) and vanish for the very central collisions which 
are cylindrical symmetric. The expansion parameters are 
v\ for directed flow and V2 for elliptic flow. 

At Bevalac energies where the directed flow in semi- 
central collisions of heavy nuclei is up to 20 degrees (lj] 
with respect to the beam line. The directed flow decrease 
with energy and is only ~ 1/20 degree at SPS energies 
At AGS and SPS energies there are now very accu- 
rate measurements for protons, positive and negatively 
charged pions as function of rapidity and transverse mo- 
mentum. Both pions and protons are found to have 
elliptic flow in the reaction plane at AGS M and SPS || 
energies, i.e., vi > for all rapidities and transverse mo- 
menta. At SPS ^2^0.1 for protons but smaller for pions 
as well as for protons at AGS. We refer to for recent 
results on flow in relativistic heavy ion collisions. 

At collision energies below a few GcV the nuclei 
shadow the collision region in the reaction plane. Conse- 
quently, the flow is "squeezed out" in the y-direction, i.e., 
V2 < 0. At ultrarelativistic energies Lorentz contraction 
of the nuclei reduce the amount of shadowing at midra- 
pidities and one finds that V2 > 0. We will in the fol- 
lowing restrict ourselves to ultrarelativistic energies and 
midrapidities where effects of shadowing are minor. 
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A. Collisionless limit/Peripheral collisions 

In peripheral collisions the nuclear overlap zone is small 
and at relativistic energies the expansion is fast. Pro- 
duced particles can therefore escape frow the collision 
zone almost without interacting with the other parti- 
cles, i.e., the system is close to free streaming and the 
collisionless limit. We can then calculate the first or- 
der correction to free streaming from particle collisions. 
This approximation is valid when the particle mean free 
paths X m fp — ((T/o)" 1 are larger than the overlap zone 
Rx,y(b). For pions and protons X m fp of order a few fm's 
for particle densities of order nuclear saturation densi- 
ties, i.e., comparable to source sizes in semi-central colli- 
sions. If we assume that particles initially are produced 
azimuthally symmetric in momentum space but not in 
coordinate space, these subsequent interactions with co- 
movers will produce an azimuthally asymmetric momen- 
tum distribution because the source is azimuthally asym- 
metric spatially. 

We shall not attempt to describe the initial hard 
nucleon-nucleon collisions, the fragmentation or particle 
production. Our starting point is some initial conditions 
at formation time tq ~ 1 fm/c and consider the sub- 
sequent scatterings between comovers described by the 
Boltzmann or Nordheim 113] equation 
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v„ • V ni = 



d 3 P2 , n da 



x [713714(1 ± 7i 1 )(l ± 7i 2 ) - 7i 1 7i 2 (l ± n 3 )(l ± 71 4 )] . (6) 

Here, n, = ?i(r i ,p i , t) is the particle distribution function 
and da/dQ, is the cms differential cross section for scat- 
tering particles 1 + 2^3 + 4. Stimulated emission and 
Pauli Blocking are including by the ± factors. 

To evaluate the effect of collisions on the distribution 
function in (|^) we assume Bjorken scaling v z = z/t ini- 
tially at invariant time tq. As usual, the space-time ra- 
pidity is 77 = (1/2) log((t + z)/(t - z)) and r = V< 2 - z 2 
is the invariant time. The transverse particle density dis- 
tribution is approximated by simple gaussians as in Eq. 
(0) . The free streaming distribution is then at later times 



n(x,p) = fo(p±,p' z ) Sj_(r ± - v_l(t - r )) . 
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Here, = (x, y) is the transverse radius, vi = (v x ,v y ) 
the transverse velocity. The local particle momentum 
distribution, /o(p), is discussed in detail in Appendix 
A; in the final result it can be rewritten in terms of final 
measured momentum distributions, dN/d 2 p t . The scaled 
longitudinal momentum 



p z = — mj_ sinh(77 — y) . 
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takes the free streaming longitudinally into account [[13| . 
Transverse expansion of the source is incorporated by 
the Galilei transformations of the transverse coordinates 



(x, y). Note, that the momentum distribution of pro- 
duced particles is azimuthally symmetric at initial pro- 
duction time To. The detailed form of /o will not be 
necessary for the evaluation of elliptic flow or HBT radii. 
It is sufficient to know the rapidity density dNj/dy of the 
scattering particles j. 

In the collisionless limit we can now insert the free 
streaming distribution function in the collision term in 
order to calculate the first order correction to the distri- 
bution function from (^). This correction provides the 
deviation from cylindrical symmetry and directly leads 
to elliptic flow, v 2 , for particle species i — Tt,p,K, as 
is evaluated in detail in Appendix A 
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Here, Vi is the particle, Vj the scatterer and vij the rel- 
ative velocity. The averages (..) refers to averaging over 
scatterer momenta pj . Since it is the momentum transfer 
in collisions, that deforms the particle momentum distri- 
bution around cylindrical symmetry, it is naturally the 
momentum transport cross section o-\ 3 r that enters. The 
elliptic flow parameter is proportional to the deformation 
and neccessarily vanishes for an azimuthally symmetric 
source. 

The other important factors in the elliptic flow are 



dNj (vijO-jj) 
dy nRxRy 
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It is the cross section times a transverse particle density 
and thus describes the effective "opacity" of the source 
for particles i scattering with particles j. 

It is remarkable that the initial time tq does not en- 
ter in (^). The ro-dependence cancels to first approxi- 
mation as explained in appendix A, due to a compen- 
sation between the densities decreasing with expansion 
time and that scatterings only lead to asymmetries as 
they pass through the source with time. The physical 
reason for the "late time" dominance in elliptic flow is 
that a particle has to travel a distance ~ R before it 
feels the source deformation, i.e., whether it moves in 
the x-direction, or in the y— direction which differs by 
a distance R y — R x ~ SR. This also explains why the 
elliptic flow is proportional to vf ± . The slow particles do 
not travel far enough to feel the deformation before the 
source has expanded - by velocity Vij t ±. 

The centrality dependence of the elliptic flow can be 
estimated from Eq. (^) . In high energy nuclear collisions 
the total multiplicity dN/dy scales approximately with 
the number of participating nucleons, Eq. (|l]) which is 
shown in Fig. 2 as function of impact parameter together 
with R x (b), R y (b) and 8(b). The resulting centrality de- 
pendence of V2 from Eq. (O), assuming that dNj/dy also 
are proportional to the number of participants, is shown 
in Fig. 2 as well. It vanishes for central and grazing col- 
lisions, has a maximum for peripheral collisions b ~ 1.5R 
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and decreases almost linearly with increasing centrality. 
A very similar behavior was found in the hydrodynamic 
calculations of Ollitrault 0. 



B. Hydrodynamic limit /Semi-central collisions 

In semi-central collisions at relativistic energies, where 
large particle densities are produced, rescatterings are 
abundant and the hydrodynamic limit (see, e.g., 0) is 
better than the collisionless — at least up to freeze-out. 

We will here attempt to extract some general features 
in the hydrodynamic limit and derive an approximate an- 
alytical formula for the elliptic flow. Instead of introduc- 
ing two independent transverse flow velocities (u x ,u y ) 
we relate the transverse flow to the spatial deforma- 
tion by the plausible assumption that the transverse flow 
is equipotential, i.e., perpendicular to and with constant 
magnitude on equi-density surfaces. This assumption is 
obvious for cylindrical symmetric sources. Also for a very 
deformed source, which appears as a slab in one trans- 
verse direction, the flow is perpendicular to most of the 
freeze-out surface. For a rectangular source the flow is 
initially also perpendicular to the surface in all directions 
but as the source expands deviations may occur at the 
corners. Concerning the constant magnitude of the trans- 
verse flow on the freeze-out surface, that may be justified 
for deformed sources by the following observation. In the 
Riemann solution to ID hydrodynamics, the flow veloc- 
ity only depends on time and the relative distance from 
the initial surface, i.e., (r — R). As shown in juj the Rie- 
mann solution is also a very good approximation to the 
cylindrical case with longitudinal Bjorken scaling with 
an additional scaling factor of (to/t) Cs for temperatures. 
For the slab or the rectangular source, the magnitude 
flow velocities will therefore only depend on the relative 
distances from the initial transverse radii, i.e., (x — R x ) 
or (y — R y ). How good this approximation is in general 
for near central relativistic heavy ion collisions should, of 
course, be checked by solving full 3D hydrodynamics. 

With the simplifying approximation of equipotential 
flow we can calculate the resulting particle spectrum for 
a locally thermal distribution e v ' u l T with longitudinal 
Bjorken flow, 



u = (7 cosh(?7), 7 sinh(?7), uj_) 



(11) 



where 7 = y/l + u\. The transverse extent is 
parametrized by (||) . With these assumptions the result- 
ing momentum distribution can be found by integrating 
over the deformed equipotential surfaces as described in 
Appendix B. The asymmetric term gives the elliptic flow 
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It is expected that the elliptic flow is proportional to the 
deformation which vanishes for central collisions. The 



relevant deformation is now at freeze-out where the radii 
which may differ from the initial transverse radii as dis- 
cussed above. 

The elliptic flow of Eqs. (|l^) depends only indirectly 
on the unknown quantities as the initial densities, equa- 
tion of state, expansion, etc. through the average trans- 
verse flow, temperature and deformation. The transverse 
flow and temperature can be extracted independently 
from transverse flow analyses by looking at the parti- 
cle mass dependence of average p± as suggested in Q 
or by measuring the apparent temperatures |16[ ], i.e., the 
inverse of the m± slopes. As shown in Appendix C the 
apparent temperature obtained from exponential fits to 
transverse mass spectra of massive particles is 



T, 



app 



T+^m(u 2 A 



(13) 



for large particle masses m and small transverse flow. Ex- 
perimentally, the apparent temperatures scales approxi- 
mately linearly with the masses of the pion, kaon, proton, 
and deuterium in central S + S and Pb + Pb collisions 
at AGS and CERN jl6| energies. From the experimental 
slopes 2dT/dm we thus obtain ~ 0.15c 2 and 0.3c 2 
for central S + S and Pb + Pb collisions at SPS energies. 

It is instructive to compare to another type of hydrody- 
namic flow parametrized by two transverse flow velocities 



u_l = {u x x/R xl u y y/R y ) . 



(14) 



The resulting particle distribution can be calculated 
along the lines of Appendix B. The resulting elliptic flow 
is for small transverse momenta 



±^((ul) - (u 2 y )) , p ± u ± ^T 

" 2= < 1(^-1)^^ , P XUX»T>- (15) 



Comparing Eq. (^5|) to the elliptic flow data for protons 
of Fig. 3 we can estimate the difference (m 2 ) — (u y ) ~ 
0.015, i.e., an order of magnitude smaller than the aver- 
age transverse flow v?,~ (u x + u y )/2 ~ 0.15 that was 
estimated from Eq. Q. The elliptic flow of Eq. (]l|) is 
not related to the spatial deformation as was the equipo- 
tential flow and therefore its predictive power is less. 



C. Comparison to elliptic flow at SPS energies 

The elliptic flow as estimated by the collisionless and 
hydrodynamic limits may be compared to NA49 data at 
SPS energies [3j taken for semicentral Pb + Pb collisions 
(b ~ R ~ 7/m). The p± dependence of V2 for pions 
and protons is measured at forward rapidities 4 < y < 5. 
We can take the free cross sections and ignore medium 
effects since the densities are small at late times when 
rescatterings lead to anisotropies (see Appendix A). The 
relevant scattering cross section is that where nucle- 
ons act as scatterers because they are heavier than e.g. 
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pions and thus deflect more. Furthermore, typical scat- 
tering cross sections cr 7rN are larger than er 71 ^ ~ lOmb 
fL7f . The former is dominated by the A-resonance. Av- 
eraging of scatterer momenta we estimate a vN ~ 30mb 



and 



..N N 



40mb. The transport cross sections are 



smaller because forward scattering without momentum 
loss should not be included. For an estimate we simply 
take the transport cross section as half of the total cross 
section. The relevant rapidity density is that of nucleons 
which we take from the NA49 experiment dNw jdy ~ 30. 
From Fig. 2 we obtain R x ~ 2/m, R y ~ 3/m and 
5 ~ 0.4 for b = R. The resulting elliptic flow is shown 
in Fig. 3 for pions and protons. The estimates in the 
collisionless limit of Eq. (||) give a reasonable description 
of both the magnitude and pj_-dependence for both pi- 
ons and protons. Furthermore, the measured elliptic flow 
is largest at midrapidities and is approximately propor- 
tional to dN/dy as predicted by Eq.(||). 

In the hydrodynamic limit transverse flow couples to 
p± and elliptic flow is therefore expected to be the same 
for all particles at the same p± . The elliptic flow can dif- 
fer if various particle species freeze-out at different tem- 
peratures or resonance decays affect the final distribu- 
tions. The elliptic flow is shown in Fig. 3 for equipoten- 



tial flow, Eq. 



with T ~ 160 MeV and 



0.15c 2 



as estimated above. An excellent fit to the magnitude 
and p± dependence of proton elliptic flow can be ob- 
tained with a deformation at freeze-out of 6 ~ 0.07. This 
value for the deformation is a factor ~ 6 smaller than the 
initial deformation of the collision zone. It is expected 
to be smaller due to expansion between initial collision 
and freeze-out but a full 3+1 dimensional hydrodynamic 
calculation is needed in order to chech the magnitude 
of the deformation as well as the validity of equipoten- 
tial flow. The deformation at freeze-out can in principle 
be extracted from HBT analyses as shown in the follow- 
ing section and would thus be an important independent 
check. 

The pion elliptic flow cannot be explained in the hy- 
drodynamic limit as it differs from proton flow and does 
not have the right p^ dependence. It is curious that 
the magnitude of the pion elliptic flow lies between the 
collisionless and hydrodynamic limit. This may indicate 
that semi-central Pb + Pb collisions at SPS energies lies 
between these two limits as one would expect since the 
particle mean free paths are a few /ermi's — compara- 
ble to the transverse source sizes for b = Rpb — 7 fm/c. 
Comparing the p^-dependence of elliptic flow for pions 
and protons would reveal the change from collisionless ex- 
pansion to hydrodynamic flow. In the collisionless limit 
the pion and proton elliptic flow differ whereas in the hy- 
drodynamic limit they should become the same as shown 
in Fig. 3. 

Calculations with the RQMD and VENUS models (see 
for detailed comparisons) give an almost rapidity in- 
dependent elliptic flow which underpredicts the SPS data 
by a factor ~ 2 — 4 at midrapidity. The smaller elliptic 
flow in RQMD at midrapidities is attributed to "preequi- 



librium softening" B. Hydrodynamic models results in 
corresponding values for vi which are ~ 4 times larger 
for an ideal (c 2 . = 1/3) pion gas and ~ 2 times larger 
when a first order phase transition at T c ~ 150MeV is 
included 0. 

It is curious that the elliptic flow is less than a few per- 
cent in relativistic nuclear collisions in comparison to the 
initial spatial deformations which are of order 5 ~ 50% in 
semi-central collisions. The reason can be understood in 
both the collisionless and hydrodynamic limits due to the 
limited time and distances that particles have to rescat- 
ter and develop asymmetric collective flow before they 
freeze-out. In other words, the scatterings that lead to 
asymmetric flow occur at late times where the source and 
scatterers have already expanded and reduced the initial 
deformation. 

The elliptic flow at RHIC and LHC energies would be 
large if one simply extrapolates from intermediate, AGS 
and SPS energies || However, at the intermediate ener- 
gies shadowing is responsible for a negative v-i (squeeze- 
out) and is still felt at AGS energies. At SPS energies 
it may also affect the forward rapidity data 4 < y < 5 
where the projectile nucleus is less Lorentz contracted 
and thus shadow. But at central rapidities at SPS and 
higher energies shadowing is minor due to the strong 
Lorentz contraction of nuclear sizes. Therefore, v-i is un- 
affected by shadowing from SPS and up in energy, i.e., 
constant rather than increase as given by a naive ex- 
trapolation from intermediate, AGS and SPS energies. 
However, since dN/dy and probably also the transverse 
flow is larger at RHIC and LHC energies, we may expect 
stronger elliptic flow. On the other hand, the expansion 
reduces the deformation and thus also the elliptic flow. 



IV. HBT 

A brief description of interferometry will be given de- 
scribing how to calculate correlation functions and HBT 
radii from a given source. The HBT radii will then be 
calculated for deformed sources - transparent as well as 
opaque. A brief description of results were presented in 
|l8| . A comparison is given and the advantages of a com- 
bined HBT and elliptic flow analysis is discussed. 

Particle interferometry was invented by Hanbury- 
Brown & Twiss (HBT) for stellar size determination Q] 
and is now employed in nuclear collisions [^0]-^4 10 Tj| . It 
is a very powerful method to determine the 3-dimcnsional 
source sizes, life-times, duration of emission, flow, etc. of 
pions, kaons, etc. at freeze-out. Since the number of 
pairs grow with the multiplicity per event squared the 
HBT method will become even better at RHIC and LHC 
colliders where the multiplicity will be even higher. 

The standard HBT method for calculating the Bose- 
Einstein correlation function from the interference of two 
identical particles is now briefly discussed. For a source 
of size R we consider two particles emitted a distance 
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~ R apart with relative momentum q = (ki — k 2 ) and 
average momentum, K = (ki +k2)/2. Typical heavy ion 
sources in nuclear collisions are of size R ~ 5 fm, so that 
interference occurs predominantly when q^h/R ~ 40 
MeV/c. Since typical particle momenta are k{^K ~ 300 
MeV/c, the interfering particles travel almost parallel, 
i.e., fci ~ k,2 — K 3> q. The correlation function due to 
Bose-Einstein interference of identical spin zero bosons 
as 7r ± 7r ± , K ± K ± , etc. from an incoherent source is (see, 
e.g., H) 



C 2 (q,K) = l ± 



/ d A x S(x,K) e 



iqx 



Jd 4 x S(x,K) 



(16) 



where S(x, K) is the source distribution function describ- 
ing the phase space density of the emitting source. 

Experimentally the correlation functions are often 
parametrized by the gaussian form 



C-2(qs,q ,qi) = 1 + \exp[-q 2 R% - q 2 R 2 

-2q q s R 2 0S - 2q q l R 2 ol ] . 



r, 2 R 2 

1i K i 



(17) 



Here, q = ki — k2 = {q s ,Qo,<li) is the relative momen- 
tum between the two particles and Ri,i — s,o,l,os,ol 
the corresponding sideward, outward, longitudinal, out- 
sideward and out-long HBT radii respectively. We have 
suppressed the K dependence. We will employ the stan- 
dard geometry, where the longitudinal direction is along 
the beam axis, the outward direction is along K_l, and 
the sideward axis is perpendicular to these. Usually, each 
pair of particles is lorentz boosted longitudinal to the sys- 
tem where their rapidity vanishes, y = 0. Their average 
momentum K is then perpendicular to the beam axis 
and is chosen as the outward direction. In this system 
the pair velocity /3j£-=K/ Ek points in the outward di- 
rection with p o = p±/m±, where m± = y/ m 2 + p\ is 
the transverse mass, and the out-longitudinal coupling 
R i vanish at midrapidity |23|]. Also R os vanishes for a 
cylindrically symmetric source or if the azimuthal angle 
of the reaction plane is not determined and therefore av- 
eraged over — as has been the case experimentally so 
far. The reduction factor A in Eq. ( |l7| ) may be due to 
long lived resonances ^,^,|), coherence effects, incor- 
rect Coulomb corrections or other effects. It is A ~ 0.5 
for pions and A ~ 0.9 for kaons. 

The Bose-Einstein correlation function can now be cal- 
culated for a deformed source. Let us first investigate 
transparent sources assuming that it is equally likely that 
a particle arrives at the detector from the front side 
of the source as from the back side, i.e., rescatterings 
through the source and opacities are ignored. As in ||, 
we parametrize the transverse and temporal extent by 
gaussians 



S(x, K) ~ S x (x, y) exp[- (T 2 ^ )2 ] e^' T 



(18) 



with longitudinal Bjorkcn flow, u = (cosh 77, 0, 0, sinhry), 
as in Eq ( pT| ) without transverse flow. Effects of trans- 
verse flow will be discussed below. The transverse radii 



R x ,R y are the gaussian radii at freeze-out, 17 is the 
freeze-out time and St the duration of emission. 

In order to calculate the correlation function of (Jlq) the 
gaussian approximations is employed (see, e.g., p2| , p3| ) 
which results in a correlation function of the form given 
in Eq. (fTlj). Inserting the source ( |l8| ) in Eq. ( |l6| ) and 
Fourier transforming we obtain the correlation function. 
The HBT radii are then obtained by comparing to the 
experimental correlation function of Eq.(17). 

It is convenient to introduce the source average of a 
quantity O defined by 



= Jd A x S(x,K)Q 
[ '~ /d 4 iS(x,K) 



(19) 



With qx ~ q • x — q-/3if t one can, by expanding Eq. ( fig 1 ) 
to second order in qiRi and compare to Eq. (0), find 
the HBT radii I2I 



Ri 

Rl 



(0* 
(0* 



Pi t) 2 ) - 

f3i t)(Xj 



Pi tf 



(20) 



P J t))-(x l -p i t)(x J -P 1 t) 



(21) 



with i,j=s^o,l. The HBT radii are a measure for the fluc- 
tuations |p5[ , variance or "length of homogeneity" |2(| 
of (xi — Pit) over the source emission function S(x,K). 
One should notice that the coordinates (x ,x s ) are ro- 
tated with respect to the (x, y) reaction plane (see Fig. 
1) by the azimuthal angle <fi between the transverse mo- 
mentum p± and the reaction plane. The beam axis is 
the same xi = z. In the local cms frame (y = 0) the pair 
velocity is P D = P± whereas P s = Pi — 0. All Rij vanish 
in this frame for cylindrical symmetric sources and for an 
azimuthally asymmetric source only R os is nonvanishing. 

For transparent sources the azimuthal dependence of 
the HBT radii has been calculated in detail by Wiede- 
mann In the longitudinal center-of-mass system of 
the pair (y = 0), the HBT radii are 



R 2 S = R 2 [1 + <5cos(20)] , 
R 2 a = R 2 [1 - <5cos(20)] + P 2 Jt 2 
R 2 os = i? 2 <Ssin(20) , 

T , 



St' 



(22) 
(23) 
(24) 

(25) 



where R 2 — (R 2 + R 2 )/2 is the average of the source 
radii squared. One notices the characteristic modula- 
tion cos(20) due to the rotation of the axes i? 0jS with 
respect to the reaction plane. Near target and projectile 
rapidities the directed flow is appreciable and leads to 
coscj) terms in Eqs. ( ^2[]23| ) [0J^]. The out- and sideward 
HBT radii show a characteristic modulation as function 
of azimuthal angle with amplitude of same magnitude 
but opposite sign. Measuring the amplitude modulation 
of Rs.o,os determines five quantities and thus overdeter- 
mines the three source parameters, namely the source 
size R, deformation 5 and and duration of emission St. 



G 



Next we consider opaque sources. In relativistic heavy 
ion collisions source sizes and densities are large and one 
would expect rescatterings. As a result particles are pre- 
dominantly emitted near the surface and arrive from the 
(front) side of the source facing towards the detector. In 
|P5| it was found that for opaque sources, where mean free 
paths are smaller than source sizes, X m f P ^R, the side- 
ward HBT radius increase whereas the outward is sig- 
nificantly reduced simply because the surface emission 
region is smaller than the whole source. As in [ p5|j27| 
Glauber absorption is introduced by adding an absorp- 
tion factor 



5 Q b s (x, K) ~ exp(— / ap(x')dx') 



(26) 



where a is the interaction cross section, p the density of 
scatterers and the integral runs along the particle trajec- 
tory from source point x to the detector. Defining the 
mean free path as \ m fp = (crp(0)) _1 , where p(0) is the 
central density, the source is opaque when X m f p <C R 
and transparent when X m f P 3> R. Calculating the corre- 
lation function for an opaque source from Eq. (16) and 
comparing to the definition of the HBT radii in Eq. (17), 
one generally obtains for small deformations, 



R 2 s =g s R 2 [l + 5cos(2^)] , 
R 2 o = 9o R 2 [l -5cos{2<t>)] + (3 2 Jt 2 
R 2 os = g os R 2 6 sm(2$), 



(27) 
(28) 
(29) 



and R 2 is unchanged from Eq. (p5|). Here g , s ,os & r e 
model dependent factors that are functions of opacity but 
independent of the deformation and can be calculated as 
for cylindrical symmetric sources. For a gaussian source 
(p oc S±), which is moderately opaque (X m f p /R — 1), a 
numerical calculation gives g s ~ 1.4 and g ~ 0.9 (see 
also p7|). For a disk source with the same rms trans- 
verse radius as the gaussian radius (i.e., Rdisk = 2i?), 
that emits like a black body (i.e., X m fp <C R), one finds 
g s = 4/3 and g Q = 4(§ - (f ) 2 ) ~ 0.2 [||. Generally, 
(.9s ~ 9o) > and the difference increases with opacity. 
Only for a completely transparent sources is g s = g . In 
all cases g os ~ g s . In Fig. 4 the HBT radii of Eqs. ( p2[ - 
|29| ) are shown for a near-central collisions (6 — 0.4) with 
a moderate duration of emission {(3 2 8t 2 /g s R 2 — 0.5) for 
various opacities X m f p /R = 0.1, 0.5, 1.0, 2.0, oo. As the 
opacity increases, g /gs decreases and therefore also the 
outward HBT radius and its amplitude. 

Comparing the HBT radii from an opaque source 
Eqs.( p7j]29| ) with those of a transparent source Eqs.(p2f- 
|24j), one notices that the amplitudes in R s and R differ 
by the amount (g s — g Q ). The modulation of the HBT 
radii with <f> provides five measurable quantities which 
over-determinates the four physical quantities describing 
the source: its size R, deformation S, opacity (g a — g s ) 
and duration of emission St, at each impact parameter. 
The azimuthal dependence of the HBT radii thus offers 
an unique way to determine the opacity of the source as 
well as the duration of emission separately. 



Experimentally, HBT analyses have not been combined 
with determination of the reaction plane yet. Conse- 
quently, the azimuthal angle <fi is averaged and the in- 
formation on three of the five measurable quantities in 
Eqs. (|27j-p9"|) is lost. From the angular averaged differ- 
ence between the out- and sideward HBT radii 



(R 2 -R 2 s )^p 2 Sr 2 ~(g s -g )R 2 



(30) 



one can only determine the sum of the positive duration 
of emission and negative opacity effect. Experimentally 
this difference is small; NA49 |l(J and NA44 Q data 
even differ on the sign. Detailed analyses of the p± depen- 
dence of the HBT radii from NA49 data within opaque 
sources |27 indicate that the sources are transparent or 
at most moderately opaque. However, the NA44 data, 
for which R ^R Sl requires an opaque source as seen from 
Eq. (|3C|). Furthermore, the p± dependence of the trans- 
verse HBT radii change if the source sizes, opacities, and 
duration of emission also are p^-dependent. 

Transverse flow may affect the out- and sideward HBT 
radii as opacity, i.e., the factors g , s may depend on both. 
For transparent sources transverse flow has been studied 
in through Eq. ( |i"l| ) assuming that the trans- 

verse flow scales with transverse distance, ~ uqt_\_/R. 
Otherwise the same transparent gaussian source as in 
Eq. ( |l8| ) was employed. To lowest order in the trans- 
verse flow both transverse HBT radii decrease by the 
same factor ~ (1 +UqTO_l/T) to leading order in uq. This 
transverse flow correction is independent of the source 
size and therefore also the deformation. Consequently, 
spatial deformations reduce the amplitudes by the same 
amount in this model for both transparent (see ||) and 
opaque sources. There are, however, box shape models 
where the transverse flow reduce R a more than R s j27| . A 
similar flow effect is found in hydrodynamic models also 
f35f although it is considerably less than the opacity ef- 
fect, i.e., when the Cooper- Frye freeze-out condition |2^] 
is replaced by the Bugaev freeze-out J2^| . Preferrably, a 
transverse flow analysis of particle spectra should be per- 
formed in order to determine the magnitude of the flow so 
that it can be separated from opacity effects. The crucial 
question, however, is whether the transverse flow is az- 
imuthally dependent, i.e., g 8>0 depend on <j> and thereby 
change the amplitudes? We can estimate this effect from 
the elliptic flow. For a simple thermal source ( Ji~8| ) with 
transverse flow as in Eq. (|i~4|), the resulting elliptic flow 
Eq. (H||) led to a very small difference for semicentral 
Pb + Pb collisions at SPS energies {u 2 } — (u 2 ) less than 
a few percent. Therefore, the azimuthally dependence 
of the flow and its effect on the amplitude in the HBT 
radii is also of that order only which is much less than 
S ~ 0.5 for semi-central collisions. Both the deformation 
and elliptic flow decrease with centrality. The conclusion 
is that besides opacity also transverse flow may affect 
the factors g 0jS - but independently of azimuthal angle. 
Therefore Eqs. (p7[-p9[) are still valid and can be used to 
extract the duration of emission unambiguously. 
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It is very important to measure the centrality or im- 
pact parameter dependence of the source sizes, deforma- 
tion, opacity, emission times and duration of emission 
in order to determine how the source change with ini- 
tial energy density. If no phase transition takes place 
one would expect that source sizes and emission times 
increase gradually with centrality whereas the deforma- 
tion decrease approximately as in Eq. (^). In peripheral 
collisions, where source sizes and densities are small and 
few rescatterings occur, the source is transparent and the 
HBT radii are given by Eqs. ( p2[|25| ). For near central 
collisions source sizes and densities are higher which leads 
to more rescatterings. Thus the source is more opaque 
and the amplitudes should differ. It would be interesting 
to observe this gradual change in the amplitudes with 
centrality. At the same time it would provide a direct 
experimental determination whether the source is trans- 
parent or opaque as well as extracting the magnitudes of 
the opacity and duration of emission separately. 



V. EFFECTS OF A PHASE TRANSITION 

If the matter created in an ultra-relativistic heavy ion 
collision undergoes a phase transition it may affect the 
elliptic flow and the modulation of the HBT radii. The 
critical energy density can be overcome either by increas- 
ing the centrality, sizes of the colliding nuclei and/or by 
increasing the collision energy as will be possible at the 
RHIC collider, y/s ~ 20 -> 200 GeV. For energy densities 
just above the critical value a mixed phase is created for 
a first order transition with zero compressibility, i.e., van- 
ishing sound speed. Even for a second order or smooth 
transition, the sound speed c 2 s = dP/de is significantly re- 
duced according to lattice calculation s |30(| . The resulting 
transverse flow is therefore reduced |15| and should lead 
to smaller elliptic flow according to Eq. (|l2|). 

Unfortunately, the elliptic flow and the modulation of 
the HBT radii both scale with the deformation and thus 
decrease linearly with centrality and vanish for near cen- 
tral collisions, where the phase transition might occur. 
However, by taking the ratio of the elliptic flow to the 
amplitude of the HBT radii v%j 8, the deformation can- 
cel and the purely geometrical effects are removed. Only 
the collective effects as transverse flow are left to first 
approximation. Since the transverse flow has been found 
to increase with the size of the collisions system at AGS 
and SPS energies, we can also expect it to increase with 
centrality in the absence of phase transitions. In the 
presense of a phase transition, however, the transverse 
flow increase very little with centrality which results in a 
"plateau" in the average transverse momentum as func- 
tion of centrality [ pl|Jl5|| . The almost constant transverse 
flow for centralities that produce a mixed phase should 
therefore according to Eq. ( |l2| ) also show up in the ratio 
vijb versus centrality such that the monotonic increase 
in transverse flow and V2 /S is replaced at the critical cen- 



trality by a plateau as function of centrality as sketched 
in Fig. 5. The plateau extends in an interval of semi- 
centrality, c\ — ► c%, where the pressure is constant due to 
a mixed phase, i.e., vanishing sound speed. 

The functional dependence of V2 /8 would be similar to 
the van Hove prediction of the temperature as function 
of rapidity density [j3l| . Also it is similar to the caloric 
curve as was recently observed in intermediate energy 
nuclear collisions |52|. The caloric curve plots the tem- 
perature versus excitation energy and displays a plateau 
as was predicted for the liquid-gas phase transition of nu- 
clear matter [ 34| . The interpretation of the caloric curve 
is, however, debated in the low energy community (see, 
e.g., H). 

If a phase transition occurs at some centrality, where 
energy densities exceed the critical value, one may also 
observe sudden changes in the physical quantities mea- 
sured in HBT analyses. The emission time and duration 
of emission increase drastically in hydrodynamic calcu- 
lations [[H leading to very large Ri and R (see also 
Fig. 5). A long lived mixed phase would also emit par- 
ticles as a black body and thus the opacity increase and 
the opacity factor g /g s decrease as shown in Fig. 5. If 
droplet formation occurs leading to rapidity fluctuations, 
one may be able to trigger on such fluctuations and find 
smaller longitudinal and sideward HBT radii p6[ . 

Alternatively, one could plot the quantities as function 
of collision energy or size of the colliding nuclei for fixed 
centrality with a similar qualitative dependence of the 
physical quantities with and without a phase transition. 
Furthermore, if an interesting change in these quanti- 
ties should occur at some centrality or collision energy, it 
would also be most interesting to look for simultaneous 
J/4" suppression [[37J, strangeness enhancement, decrease 
in directed, elliptic or transverse flow, or other signals 
from a phase transition. 



VI. SUMMARY 

In summary, measuring the reaction plane in relativis- 
tic nuclear collisions, the elliptic flow and HBT radii si- 
multaneously gives a detailed description of the source 
at freeze-out as well as how the source expands in space 
and time from the initial collision and up to freeze-out. 

The elliptic flow was calculated in the collisionless limit 
from the Boltzmann equation with a collision term as 
well as in the hydrodynamic limit assuming equipoten- 
tial flow. The collisionless limit is expected to be valid for 
peripheral collision where few rescatterings are expected. 
Alternatively, the hydrodynamic limit is expected to be 
valid for central collision where many rescatterings are 
expected. For small transverse source deformations the 
elliptic flow is proportional to the deformation in both 
limits. Detailed cascade and 3+1 dimensional hydrody- 
namic calculations are clearly needed of both the radial 
and elliptic flow as well as HBT radii for non-central nu- 
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clear collisions with and without phase transitions. 

The estimates in the collisionless limit of Eq.(^) give 
a reasonable description of both the magnitude and p±- 
dependence of the elliptic flow measured at SPS energies 
in Pb + Pb semicentral collisions for both pions and pro- 
tons. Furthermore, the measured elliptic flow is largest at 
midrapidies and is approximately proportional to dN/dy 
as also predicted by Eq.(||). In the hydrodynamic limit 
transverse flow couples to p± and elliptic flow is therefore 
expected to be the same for all particles at the same p± . 
For equipotential flow given by Eq. ( |l2| ) can also describe 
the proton elliptic flow , however, with a deformation is 
a factor ~ 7 smaller than the initial deformation of the 
collision zone, ft is expected to be somewhat smaller due 
to expansion between initial collision and freeze-out. The 
measured elliptic flow for pions lies between the collision- 
less and hydrodynamic limit as could be expected since 
the particle mean free paths are a few fermi's — com- 
parable to the transverse source sizes for b = Rpb — 7 
fm/c. Comparing the pj_-dependence of elliptic flow for 
pions and protons reveals distinctively the change from 
collisionless expansion to hydrodynamic flow. In the col- 
lisionless limit they differ (see Fig. 3) whereas in the 
hydrodynamic limit they should become the same. 

The modulation of the HBT radii with azimuthal an- 
gle between the reaction plane and particle transverse 
momenta can be exploited to obtain source sizes, defor- 
mations, life-times, duration of emission and opacities 
separately. HBT radii provides important space-time in- 
formation which complements the momentum space in- 
formation from particle spectra on elliptic flow during 
expansion and freeze-out. The HBT radii can distin- 
guish between opaque sources Eqs.( p7H29| ) from transpar- 
ent source Eqs.( p2| - p4| ) because the amplitudes in R s and 
R Q differ. The modulation of the HBT radii with (f> pro- 
vides five measurable quantities which over-determinates 
the four physical quantities describing the source: its size 
R, deformation 8, opacity (g a — g s ) and duration of emis- 
sion St, at each impact parameter. The azimuthal de- 
pendence of the HBT radii thus offers an unique way to 
determine the opacity of the source as well as the dura- 
tion of emission separately. 

Tracking these physical quantities with centrality will 
provide detailed information about the source created in 
relativistic nuclear collisions and may reveal a phase tran- 
sition as sketched in Fig. (5). 
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APPENDIX A: Evaluation of the Collision In- 
tegral 

In order to calculate the elliptic flow, we need the par- 
ticle distribution integrated over space for normalization 
With the free streaming distribution function 
we find at any time t 




dyd 2 p\ 



= (2vr 



,-3 



m± cosh(?7 — y)n(r, p, t)rdr]d 2 r± 



(2tt) d r / dp z f (p_L,p z ). 



(31) 



The change of particle momenta due to collisions is 
the change in the particle distribution function integrated 
over space and time 

dANi Ei [ (dn 
~dt 



dt d 3 r . 



(32) 



dyd 2 p± (2tt) 3 j coll 

In the collisionless limit we can now insert the free 
streaming distribution function (Q) in the collision term 
in order to calculate the first order correction to the dis- 
tribution function. As will be shown below, the system 
expands rapidly to low densities, where collisions mainly 
contribute, and therefore the final state factors (1 ± ni) 
from stimulated emission/Pauli blocking can be ignored. 
For small angle scatterings the loss and gain terms in 
(n3?i4 — nxn-i) nearly cancel leading to a suppression fac- 
tor at forward and backward angles ~ (1 — cos 2 9). Fur- 
thermore, little momentum is lost in small angle deflec- 
tions leading to little deformation in momentum space. 
We therefore replace the cross section by an angular and 
energy averaged transport cross section cr tr and keep only 
the loss term oc n\U2 

dANi E x f , 2 , , f d 3 p 2 



dyd 2 pj 



(2tt) e 



d r±TdTdij 



(2tt) : 



■ v 12 <Jtr nin 2 . 



(33) 

where we have changed variables to the space-time ra- 
pidity and invariant time. 

At late times we can utilize that the free streaming 
distribution longitudinally can be approximated by 

fo(px,P z ) - ^rrf^Siy-rj) J dp z f (p±,p z ) . (34) 

where the latter integral can be eliminated by use of Eq. 
(|3f|). The resulting free streaming distribution function 
of Eq. becomes 



n{x,p) ~ S ± (r ± - v ± r ) 8{y - <q)- 



/ m± dyd z p± 

First we deal with the integration over transverse co- 
ordinates 







= J dxdy Sj_(r± - Vi±t') S±(r± - v 2±T ') 



4:TrR x R y 



exp 



{Vlx - V 2x f - — (Vly - V 2y f 



with t' = t — tq. The transverse particle momentum 
or velocity with respect to the reaction plane (x-axis) is 
vu_ = (vi x ,viy) — «i_l(cos <f>, sin 0), in terms of the az- 
imuthal angle 0. Assuming deformations are small we 
obtain by expanding in S 



u = 



1 



4irR r R 



■ exp 



v 

4i? 2 



'2 ,,2 



T V 



12 



4i? 2 _ 

5cos(20) + 0{8 2 ) 



(36) 



16irR x R., 



dy vL 



(39) 



APPENDIX B: Evaluation of equipotential flow 
in semi-central collisions 

In this appendix we evaluate the elliptic flow for 
an asymmetric source with equipotential flow which is 
parametrized spatially as a gaussian in transverse direc- 
tions. We assume particles are thermally distributed with 
local flow velocity as in Eq. ([ll]) . By assumption the 
transverse flow is constant on the freeze-out surface which 
is determined by 



where v± 2 — |vu_ — v 2 j_| is the relative velocity. At 
this point we make the important observation that the 
asymmetric term proportional to cos(20) is weighted by 
a factor r' 2 = (r — to) 2 , i.e., scatterings at early times 
contribute little to azimuthal asymmetries in the mo- 
mentum distribution whereas they are important around 
r ~ 2R/V-L2- This justifies Eqs. © and © as well 
as the neglect of the Bose and Fermi factors (1 ± m) in 
the collision term as densities are low at late times. The 
physical reason for the "late time" dominance in elliptic 
flow is that particles have to travel a distance ~ R before 
it feels the source deformation, i.e., whether it moves in 
the x— direction, where R x < R, or in the y— direction, 
where R y > R. 

Since tq <C R/v±, the time integral becomes 



dr 5 vi I . , s 

— 1± = — - — — — ^— cos(20) + constant , (37) 

t 8nR x R y vf 2 



where the constant is independent of azimuthal angle <f>. 
Finally, we carry out the p 2 and dn integrals in Eq. j3^ ) 
using ( p7\ ) with the result 

^AAx = ^ dN 2 dN! S v 2 1± 
dyd 2 p± r dy dyd 2 p± 8ttR x R v (vf 2 ) 

+ constant , (38) 

where the averaging over scatterer momenta is indicated 
by (..). Originally we only included loss terms in the 
Boltzmann equation and thereby replaced the cross sec- 
tion by the transport cross section. The gain terms will 
not affect the asymmetric term but will cancel the con- 
stant term in Eq. ( |38| ) when the number of particles is 
conserved. If there is net particle absorption or produc- 
tion, v 2 should be multiplied by the ratio of initial and 
final particle number. 

We can now generalize to several kind of scatterers by 
replacing 1 by a given particle species i as pions, protons, 
etc. and replace the scatterer 2 by a sum of scatterers 
j. By dividing Eq. ( p8| ) by dNi/dyd 2 p± we obtain by 
comparing to the definition, Eq. (||), of elliptic flow 



x2_ y^ 

Rl + R 2 y 



(40) 



i.e., u j_ is a function of a but independent of the az- 
imuthal angle 0' = tan _1 (y/x). Also, the transverse flow 
velocity is perpendicular to the elliptic surface 



= u±n = uj_(a) 



(i? 2 cos 0',i? 2 sin (j)') 
R* cos 2 0' + R x sin 2 



(41) 



In calculating the azimuthal dependence only the 
transverse directions are important as normalizations 
cancel in Eq. (||). Thus 



dN 



J adaS±(a) J d(j>' exp(uj_ • p±/T) . 



dyd?pi 

Expanding for small deformations we obtain 



dN 



dyd 2 p ± 



<x (M— [l 



5cos(20)])) 



(42) 



(43) 



where the average (..) refers to radial average over a of 
the transverse flow. Iq is the Bessel function of imaginary 
argument which has the limits 



Io(x) 




(44) 



Expanding for small deformations and comparing to the 
definition of elliptic flow (|J) we finally obtain 



V-2 




1)5 



p±u± -C T 
p±u± >• T 



(45) 



APPENDIX C: Apparent Temperatures and 
Transverse Flow 
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Transverse flow affects the measured p± and m± 
slopes. The apparent temperatures, obtained by fitting 
the particle spectra by dN/d 2 p± oc exp(— m±/T app ), are 
larger than the intrinsic ones, ft is difficult to deter- 
mine the intrinsic temperature and the transverse flow 
separately from p± slopes of pions alone |3q| . Recent 
measurements Jl6| of apparent temperatures for various 
massive particles, n, K,p, d, 3 He, etc, may allow us to 
estimate the transverse flow uniquely as will now be de- 
scribed. 

Assume a thermal source in two dimensions with in- 
trinsic temperature T at freeze-out and (transverse) flow 
Uj_ locally, i.e., n ~ exp{p ■ u/T). For small transverse 
deformations the flow is almost azimuthally symmetric 
and its average value is equal to the cylindrical symmet- 
ric value, u±(r). The distribution of particles is thus 



n(r,p, t) ~ exp[ 



-7m j 



U_l 



T 



(46) 



In three dimensions the source is further complicated by 
the thermal factor exp(— mj_ cosh(y — rj)/T). After in- 
tegrating over longitudinal direction or rj, however, the 
result are very similar when m± 3> T. 

The dependence on the angle 6 between pj_ and in 
Eq. (46) is crucial. Integrating over transverse coordi- 



nates gives the distribution 

dN m -i-M 
-p— cx exp[-7— \{Ia{pj_u_L/T)) , 



(47) 



where (...) refers to radial average of the transverse flow 
and Iq the Bessel function (see Eq.(|44|). 

The apparent temperature defined as the inverse m± 
slope becomes 



T a 



dm_ 
T + 



>n(^)]-' 



*d 2 p± 



Pa_uj_ < T , 
p±u ± > T , 



m±Uj_ <C T 
p±^> m 



(48) 



At small p± the result is the expected one when the ki- 
netic energy of flow is added to the thermal energies in 
two dimensions. At large p± and flow one instead ob- 
tains the blueshift formula J39|. Experimentally, the ap- 
parent temperature is determined by exponential fits to 
dN/d 2 p± in a certain region of p± and it may therefore 
differ somewhat from (fi"8|). 

One should notice that in experimental fits to parti- 
cle spectra the apparent temperatures are parametrized 



as T, 



app 



T + m/3 . Consequently, this flow parameter 



differs from ours by a factor: 1 — (u 2 l }/2. 
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FIG. 2. Transverse radii of nuclear overlap, deformation, number of participants and elliptic flow parameter (see Eq.(J9j) and 
text) versus impact parameter. 
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FIG. 3. Elliptic flow V2 for pions and protons versus transverse momentum. Data from NA49 |3j. Curves represent hydrody- 
namic limit of Eq.([t2]) (full curve) and collisionless limit of Eq.(|^) for pions (dashed) and protons (dotted) (see text for details). 



2.0 



1.5 



CC 

if) 

CO _ 

CC 



0.5 



— RosW 

— R 2 /g s R 2 




1.0 



0.0 ^ 
0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



§1% 



FIG. 4. HBT radii vs. angle between reaction plane and transverse particle momenta. The HBT radii are normalized to the 
angle averaged sideward HBT radius squared, g s R 2 , for deformed 5 — 0.4 source with duration of emission 0i5r 2 / g s FI? — 0.5 
and with various opacities. The sideward HBT radius (full curve) is then the same for both transparent (Eq. (E2h) and opaque 
(Eq. (^)) sources and likewise for the out-side HBT radii (Eqs . (|i^ ) and (^), dashed curve). The outward HBT radii 
are shown with chain-dashed curves for a gaussian source (Eq. ([ii^and (^) with various opacities (from below and up): 
\ mfp /R = 0.1, 0.5, 1.0, 2.0, oo. 
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FIG. 5. Skematic behavior of the centrality, dN/dy or Et dependence with and without a phase transition. From top to 
bottom the curves show the ratio of elliptic flow to deformation V2/6, opacity g /g s , and duration of emission St. The expected 
behavior without (full curves) and with (dashed curves) a first order phase transition occurring in a region of semi-centrality 
from cl to c2 is described in the text. 
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